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The electronic, magnetic and transport properties of Fe intercalated 2H-TaS2 have been investi¬ 
gated by means of the Korringa-Kohn-Rostoker (KKR) method. The non-stoichiometry and disor¬ 
der in the system has been accounted for using the Coherent Potential Approximation (CPA) alloy 
theory. A pronounced influence of disorder on the spin magnetic moment has been found for the 
ferro-magnetically ordered material. The same applies for the spin-orbit induced orbital magnetic 
moment and magneto-crystalline anisotropy energy. The temperature-dependence of the resistivity 
of disordered 2H-Feo.28TaS2 investigated on the basis of the Kubo-Stfeda formalism in combination 
with the alloy analogy model has been found in very satisfying agreement with experimental data. 
This also holds for the temperature dependent anomalous Hall resistivity pxy{T). The role of ther¬ 
mally induced lattice vibrations and spin fluctuations for the transport properties is discussed in 
detail. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

Transition metal dichalcogenides are very attractive 
materials both from a fundamental point of view as well 
as concerning potential technological applications. They 
are formed by well separated trilayers of the type MX 2 
(M: transition metal, X: S, Se, Te), that determine their 
quasi-2D structure leading to various interesting physical 
properties. These are for instance strongly anisotropic 
transport properties [1] or charge-density-wave (CDW) 
instabilities leading to the formation of CDW phases 
which can coexist at low temperature with supercon¬ 
ductivity [2-4]. Many discussions in the literature con¬ 
cern the opto-electronic properties of the semiconducting 
dichalcogenides and their relation to the features of the 
energy gap as width and direct or indirect character [5, 6]. 

Most of the transition metal dichalcogenides are non¬ 
magnetic. However, they allow the intercalation of mag¬ 
netic atoms or molecules between the X-M-X trilayers. 
This leads to an interesting class of magnetic materials 
having quasi-2D properties [7-10] which can be varied 
via the amount and type of intercalated atoms. This is 
demonstrated in the present theoretical work on Fe inter¬ 
calated 2H-TaS2, where the Fe atoms occupy the octahe¬ 
dral holes between prismatic MX 2 trilayers (see Fig. 1). 
In fact, Fe intercalated 2H-TaS2 can be seen as a proto¬ 
type material that has been investigated experimentally 


* Sergiy. Mankovskyy@cup. uni-muenchen. de 


already some years ago as the closely related intercalated 
2H-structured systems NbSe 2 , TaSe 2 , NbS 2 , and TaS 2 
with the magnetic Sd-metals Mn, Fe, Co and Ni as an 
intercalate [11-14]. As was shown, a common feature of 
these materials is the trend to create two-dimensional or¬ 
dered 2x2 and x structures within the Sd-layer 
for the concentration a; = 1/4 and 1/3, respectively, of 
the intercalate. 

Depending on concentration and intercalation atom 
type, these materials exhibit for the ordered phase fer¬ 
romagnetic (FM) [8] or anti-ferromagnetic (AFM) [11] 
magnetic order at low temperatures. The AFM order 
has been observed for the intercalated NbS 2 systems, 
e.g., Fei/ 3 NbS 2 , Coi/ 3 NbS 2 , and Nii/ 3 NbS 2 , while in the 
case of TaS 2 based alloys intercalated by Sd-metals FM 
order has been found [8] for Cri/ 3 TaS 2 (Jc = 116 K), 
Mni/ 3 TaS 2 (Tc = 70 K), and Fei/ 3 TaS 2 (Tc = 35 K). 
In the case of non-stoichiometric, disordered systems the 
Curie temperature depends in addition on the concentra¬ 
tion of the 3d element [9, 10]. A strong dependence of 
Tc on the Fe concentration was found in particular for 
Fe,rTaS 2 [9] with Tc = 90, 163, 90, and 55 K for a; = 0.2, 
0.26, 0.28, 0.34, respectively. At x = 0.45 the system 
becomes AFM with a Neel temperature Tn = 85 K [10]. 

Ordered Fei/ 4 TaS 2 possesses a very pronounced out- 
of-plane magneto-crystalline anisotropy (MCA), i.e. the 
easy axis is along the c-axis. The extremely high 
anisotropy field Ba of about 60 T [15, 16], that is by an 
order of magnitude higher than observed for Mni/ 4 TaS 2 
and Mn 4 / 4 NbS 2 [7] (about 0.5 T), leads to an Ising-type 
behavior of the Fe magnetic moment. 
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The transport properties of magnetic intercalated tran¬ 
sition metal dichalcogenides have been investigated by 
various authors. Examples for this are the electrical re¬ 
sistivity and anomalous Hall effect (AHE) [11] studied 
for alloys based on 2H-NbS2. The magneto-transport and 
superconducting properties of disordered Ee intercalation 
compounds based on NbSe 2 , TaSe 2 and TaS 2 have been 
investigated by Whitney et al. [17] for the dilute regime 
{x < 0.1). Recently, the results of magneto-transport 
measurements for Fei/ 4 TaS 2 have been reported by Mo- 
rosan et al. [18]. This material shows a strong anisotropy 
concerning the magnetization as well as the resistivity. 
Checkelsky et al. [15] focused on the temperature de¬ 
pendence of the resistivity of Fei/ 4 TaS 2 . These authors 
showed that the characteristics of the AHE for T > 50 K 
cannot be explained by the Karplus-Luttinger or Berry- 
phase mechanism. Instead it was concluded that it is 
governed by scattering processes connected with temper¬ 
ature induced phonons and magnons. These authors also 
discussed the unconventional behavior of the magneto¬ 
resistance as a function of temperature. 

A very large magneto-resistance (MR) has been found 
recently for disordered Feo. 28 TaS 2 single crystal by Hardy 
et al. [19], which is nearly 100 times stronger than that 
observed for the ordered Fei/ 4 TaS 2 compound. The au¬ 
thors point out the crucial role of spin disorder and strong 
spin-orbit coupling in the system, that can be exploited 
to create materials with large MR. 

In the present work we will focus on TaS 2 intercalated 
with Fe as a prototype material for magnetic intercalated 
transition metal dichalcogenides. While several experi¬ 
mental studies on the magnetic and transport properties 
of this system can be found in the literature, only few 
theoretical investigations [20, 21] have been done so far 
with the focus on its electronic structure assuming an 
ordered state. 


II. COMPUTATIONAL DETAILS 


lated 2H-TaS2 as shown in Fig. 1 have been taken from 
the experiment [10, 18]. For the angular momentum 



FIG. 1. The structure of the investigated Fe intercalated 2H- 
TaS 2 system showing the position occupied by the Fe atoms 
and empty spheres (ES) in between the S-Ta-S trilayers ac¬ 
cording to the Fe concentration. 

expansion of the Green function a cutoff of Imax = 3 
was applied. This is in particular important for the 
transport calculations that were based on the Kubo- 
Stfeda formalism [28] that among others gives access 
to the AHE. [29] To deal with the temperature depen¬ 
dence of the transport properties a scheme has been used 
that is based on the alloy analogy model and that ac¬ 
counts for thermal lattice vibrations [30] as well as spin 
fluctuations. [31] Within this approach the spin fluctua¬ 
tions are represented by a temperature-dependent quasi¬ 
static spin configuration corresponding to the adiabatic 
approximation. [32] The spin configuration used as an in¬ 
put for the transport calculations may be determined 
for example by performing Monte Carlo simulations or 
deduced from the experimental temperature-dependent 
magnetization. [31] Here, the later scheme has been ap¬ 
plied using experimental data from Ref. 19. 


The present theoretical investigations on the magnetic 
and transport properties of Fe intercalated 2H-TaS2 are 
based on first-principles electronic structure calculations 
which have been performed using the fully relativistic 
KKR Green function method.[22, 23] The combination 
with the coherent potential approximation (CPA) alloy 
theory allowed to deal with disorder on the Fe sub-lattice 
as well as the impact of finite temperatures on transport 
properties on the basis of the alloy analogy model (see be¬ 
low). The self-consistent calculations have been done in 
the framework of the local approximation to spin density 
functional theory (LSDA) using the parametrization for 
the exchange and correlation potential as given by Vosko 
et al. [24]. Correlation effects going beyond the level of 
LSDA have been accounted for by means of the LDA-I-U 
scheme [25-27] using for Fe the parameters U = 4.5 eV 
and J = 0.7 eV [21] throughout. 

The parameters specifying the structure of Fe interca¬ 


III. RESULTS 
A. Electronic structure 

At low temperatures, pure 2H-TaS2 exhibits a charge 
density wave (CDW) instability [33] driven by the Fermi 
surface nesting mechanism, bringing the system into the 
CDW state below Tcdw = 80 K [3]. Intercalation with 
Fe obviously strongly modifies the Fermi surface due to 
Fe-related electronic states resulting in a shift of the 
Fermi energy and a smearing of the electronic bands due 
to the disorder within the Fe layers. This suppresses the 
CDW instability and the intercalated system shows con¬ 
ventional metallic behavior. 

These findings are reflected by the results of calcula¬ 
tions on the electronic structure. Fig. 2 shows the corre¬ 
sponding spin-integrated Bloch spectral function (BSF) 
A{k,E) calculated for 2H-TaS2, disordered Feo. 25 TaS 2 , 
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as well as ordered Fei/ 4 TaS 2 . To avoid the use of an ex- 



FIG. 2. Spin-integrated Bloch spectral function A{k, E) along 
high symmetry directions of the Brillouin zone calculated 
for: (a) 2H-TaS2, within LSDA; (b) ferromagnetic disordered 
Feo. 25 TaS 2 , within LSDA-I-U; and (c) ferromagnetic ordered 
Fei/ 4 TaS 2 , within LSDA-fU. In the case of ordered com¬ 
pounds, TaS 2 and ordered Fei/ 4 TaS 2 , the BSF was calculated 
for a small imaginary part Q'(-E) = 0.001 Ry of the energy E. 

tremely fine mesh of k- and ill-points these calculations 
were done using a small imaginary part of the energy E in 
the case of the ordered compounds TaS 2 and Fei/ 4 TaS 2 . 
The resulting BSF for TaS 2 corresponds essentially to 
the dispersion relation E{k) given by Blaha [34] with 
the differences primarily to be ascribed to the impact of 
the spin-orbit coupling that is accounted for within the 
present fully relativistic calculations. Intercalation of Fe 
leads for the disordered case, apart from the exchange 
splitting, primarily to a broadening of the bands as can 
be seen from Fig. 2 b). For ordered Fei/ 4 TaS 2 , on the 
other hand, the band structure gets much more complex 
due to the increase of the size of the unit by a factor of 
4 and the occurrence of new Fe-related states (see Fig. 2 
c)). 

The various features in the BSF in Fig. 2 are also re¬ 
flected by the component and spin-resolved density of 


states (DOS) shown in Fig. 3 for Fe (a), Ta (b) and S 
(c) in pure TaS 2 , disordered Feo. 25 TaS 2 as well as or¬ 
dered Fei/ 4 TaS 2 . The rather broad Fe energy bands 



FIG. 3. Fe (a), Ta (b) and S (c) DOS in pure TaS 2 (within 
LSDA) and ferromagnetic disordered Feo. 25 TaS 2 as well as or¬ 
dered Fei/ 4 TaS 2 (within LDA-I-U). In the cases of Ta (b) and 
S (c) the thin solid and dashed red lines represent the DOS 
for two in-equivalent sites, a and b, in ordered Fei/ 4 TaS 2 . 

(Fig. 3 (a)) indicate a significant hybridization with the 
electronic states of the host atoms, especially with the 
neighboring S atoms. This can also be concluded from 
the partial DOS of S, shown in Fig. 3 (c). Obviously, 
disorder in the Fe sub-lattice leads for Fe itself to quite 
pronounced changes when compared to the ordered case; 
apart from a broadening and resulting smearing of the 
spectral features in the Fe DOS one notices a substantial 
redistribution of spectral weight together with an appar¬ 
ent change in the exchange splitting. For the S and Ta 
sublattices, on the other side, disorder leads first of all to 
a smearing of the DOS curve. Interestingly the impact 
of disorder for Ta is at least as pronounced as for the S, 
being next-nearest neighbor to Fe (see Fig. 1). 
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B. Magnetic moments 

Calculations of the magnetic moments for ferromag¬ 
netic disordered 2 H-Fea;TaS 2 have been done first using 
the GGA approach for the treatment of exchange and 
correlation. As this led to results that were too low 
when compared to experiment (see below) the LDA-I-U 
method was used instead with the corresponding param¬ 
eters given in section 11. As can be seen from Fig. 4 
this led for disordered 2 H-Fea;TaS 2 to a spin magnetic 
moment that increases nearly monotonously from about 
2.42 to 3.05 fj,B when the Fe concentration x increases 
from 0.05 to 0.5. On the other hand, the relatively large 



FIG. 4. Spin (circles) and orbital (squares) magnetic moments 
of Fe in ferromagnetic 2H-Fea:TaS2 as a function of the Fe 
concentration. Results for the disordered state are given by 
full symbols while those for ordered 2 H-Fei/ 4 TaS 2 are given 
by open symbols. 

spin-orbit induced orbital moment decreases from about 
0.82 to 0.29 fj,B- In the case of ordered 2 H-Fei/ 4 TaS 2 
quite pronounced changes compared to the disordered 
state are found; upon ordering the spin moment increases 
from 2.76 to 3.64 /tb while the orbital moment decreases 
from 0.48 to 0.16 /tb- For the GGA calculations the corre¬ 
sponding changes are 1.93 to 3.1 fiB for the spin moment 
and 0.12 to 0.11 /iB for the orbital moment, respectively. 
The increase of the spin magnetic with ordering is quite 
common and can often be associated with the DOS of the 
ordered state that is more structured and less broad than 
for the disordered state (see Fig. 3). For the same reason 
one has usually also an increase of the orbital magnetic 
moment. On the other hand, the orbital magnetic mo¬ 
ment depends much more on the details of the electronic 
structure in the vicinity of the Fermi energy than the 
spin magnetic moment. This together with the relatively 
small width of the Fe sub-band seems to be the reason for 
the observed pronounced decrease for the orbital moment 
upon ordering. 

The LDA-bU-based results differ quite substantially 
from that of Ko et al. [21] who used the GGA-I-U scheme 
with the values for U and J as given in section II and 
spin-orbit coupling treated as a perturbation. This ap¬ 
proach led for the spin and orbital magnetic moment to 


the values 2.95 and 1.0 /tb, respectively. These results 
imply that in particular the orbital magnetic moment 
depend quite sensitively on the treatment of exchange 
and correlation as well as spin orbit coupling. 

Experimental work on ordered 2 H-Fei/ 4 TaS 2 led to a 
total magnetic moment per Fe atom of 3.9 /tb and a Curie 
temperature of Tq = 160 K.[18] This is in close agree¬ 
ment with more recent work by Checkelsky et al. [15] 
(Tc = 160 K) as well as Choi et al.[16]. The later au¬ 
thors did measurements on ordered samples obtained as 
grown (AG) and by slow cooling (SC) leading to a total 
moment of 4.0 fiB and a Curie temperature of Tc = 156 
and 159 K, respectively. Additional measurements on 
quenched (Q) disordered samples led to a substantially 
lower Curie temperature of Tq = 104 K. 

On the basis of XMCD measurements at the L 2 , 3 -edges 
of Fe the ratio of the spin (rngpin) to the orbital (worb) 
moment has been estimated to be TOorb/wspin = 0.33 
[21]. Assuming an ionic configuration for Fe with a spin 
moment of 4.0 /iB an orbital moment of 1.33 /tb was sug¬ 
gested. Using instead the value 4.0 /tB for the total mo¬ 
ment as given above an orbital moment of 1.0 /(b results 
from the analysis of the XMCD measurements. 

When comparing the experimental results with the 
theoretical ones given in Fig. 4 one finds that the spin 
magnetic moment calculated for the 2 H-Fei/ 4 TaS 2 com¬ 
pound is in a good agreement with the experimental data 
while the orbital one seems to be too low. This might be 
to a large extent due to the pronounced dependence of 
the orbital magnetic moment on the computational de¬ 
tails that was discussed above. 


C. Magneto-crystalline anisotropy 

The magneto-crystalline anisotropy energy AUmca 
defined as the difference in energy for the in-plane and 
out-of-plane orientation of the magnetization has been 
calculated by means of the so-called torque method [35]. 
Ai?MCA(ic) determined this way for disordered Fea;TaS 2 
as a function of the Fe concentration x and expressed 
with respect to the formula unit (see below) increases 
monotonously with x. To compare with experimen¬ 
tal data Ai?MCA(a^) was converted into an anisotropy 
field Ba = AUMCA/ATf.u., where the magnetic mo¬ 
ment per formula unit (Mf.u.) is given by that per Fe 
atom via Mf.u. = xMpe- Using this representation 
for the magneto-crystalline anisotropy Ba{x) decreases 
monotonously with x as can be seen in Fig. 5. For 
X = 0.25 calculations have been done in addition assum¬ 
ing an ordered compound Fei/ 4 TaS 2 . As often found, 
ordering increases the anisotropy field in a quite pro¬ 
nounced way. Astonishingly, this consequence of ordering 
for A£'MCA(a^) is reversed when compared to the situa¬ 
tion for the orbital magnetic moment. Comparison of 
the theoretical result for Ba with the experimental value 
for a; = 1/4 confirms the large anisotropy field found in 
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FIG. 5. Magneto-crystalline anisotropy in the Fe-intercalated 
TaS 2 vs Fe concentration. The diamond represents the MCA 
energy measured at 2 K [15, 16] (anisotropy field Ba ~ 60 T) 


experiment. Nevertheless, the theoretical value for the 
ordered as well as for the disordered state are well above 
the experimental value. This indicates again that a cer¬ 
tain amount of disorder may be present in the sample 
investigated. 


D. Electrical resistivity 

The temperature-dependent longitudinal electrical re¬ 
sistivities, Pxx{T) = yOyy(T) and pzz(T), calculated for 
ferromagnetic disordered 2 H-Feo. 28 TaS 2 using the scheme 
described in section II are shown in Fig. 6 (a). For the in¬ 
plane longitudinal electrical resistivity Pxx(T') results are 
given in addition that account only for the lattice vibra¬ 
tions (vib) and spin fluctuations (fluct), respectively. For 
T = 0 K all individual p^^{T) curves represent the resid¬ 
ual resistivity due to the random disorder in the Fe sub¬ 
lattice. Keeping the magnetic spin configuration fixed 
and dealing only with the impact of thermal lattice vi¬ 
brations on the resistivity a pronounced and monotonous 
increase with temperature is observed for the correspond¬ 
ing curve p^^{T) (up triangles). Considering on the 
other hand (down triangles) that accounts for 

the spin fluctuations only one finds a rather weak increase 
with temperature. Accounting for both temperature de¬ 
pendent scattering mechanisms simultaneously one finds 
for the resulting curve p^^{T) (full squares) only minor 
changes compared to p^{T). As found before for fer¬ 
romagnetic bcc-Fe [31], this finding clearly shows that 
and p^^{T) are not additive, i.e. the Matthiesen 
rule is violated. Comparing the calculated p^y^(T) curve 
with the available experimental data [9, 20] one can see 
that the experimental and theoretical residual resistivi¬ 
ties agree rather well. The fact that the theoretical values 
are somewhat higher may indicate that the assumption 
of fully random disorder on the Fe sub-lattice for the cal¬ 
culations is not fully justified; i.e. there might be some 
order present in the experimentally investigated samples. 

The various curves in Fig. 6 (a) clearly show that 



FIG. 6. (a): Temperature-dependent longitudinal resistivi¬ 

ties, Pxx(r) = Pjy{T) and pzz{T) (filled and open squares) 
for ferromagnetic disordered 2H-Feo.28TaS2. The dots [9] 
and the dashed line [20] represent corresponding experimen¬ 
tal data, (b) Temperature-dependent transverse resistivity, 
Pxy(r). Results for pxx(T) and pxy(r) that were obtained ac¬ 
counting only for lattice vibrations and spin fluctuations are 
represented in both cases by up and down triangles, respec¬ 
tively. 


the thermal lattice vibrations are primarily responsible 
for the temperature dependence of Pxx(^) including the 
change of its slope close to the Curie temperature. For 
that reason this peculiar behavior cannot be ascribed to 
the temperature induced magnetic disorder as it was as¬ 
sumed before. [9, 20] Obviously, the temperature depen¬ 
dence of Pxx(^) primarily to be ascribed to the in¬ 
creasing smearing of the energy bands with temperature 
due to electron scattering caused by lattice vibrations. 
For low temperature, this results in a rather fast increase 
of the resistivity due to a corresponding increasing cross 
section for the inter-band scattering. For higher temper¬ 
atures, a saturation of the number of channels for the 
inter-band scattering seems to occur leading finally to a 
rather weak increase of PxxiT) with temperature. 

As can be seen in Fig. 6 (a), the calculated out-of-plane 
resistivity Pzz{T) for T = 0 K is larger than its in-plane 
counterpart Pxx{T) reflecting the 2D-character of the sys¬ 
tem. The finding that this difference is relatively weak 
is due to the disorder in the Fe sub-lattice that causes 
in both cases a large residual resistivity. With increas- 
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ing temperature this difference diminishes monotonously 
being nearly absent for the highest temperature consid¬ 
ered here (T = 125 K). This behavior of the resistivity 
has been observed and discussed before for example for 
the compounds NbaSn and NbsSb.[36-38] The fact that 
the p 2 z(T)-curve is not shifted rigidly against that for 
Pxx(T') clearly shows once more that the contributions to 
the resistivity due to chemical and thermal disorder are 
not simply additive. This non-additive behavior together 
with the fact that the temperature dependence for Pxx{T) 
as well as for PzziT) is obviously dominated by the ther¬ 
mal lattice vibrations leads obviously to the decreasing 
anisotropy of the resistivity with increasing temperature 
seen in Fig. 6 (a). 

The temperature-dependent transverse electrical resis¬ 
tivity Pxy(T') for disordered 2H-Fea;TaS2 is shown in Fig. 
6 (b). For T = 0 K the various theoretical curves are 
determined only by the chemical disorder on the Fe sub¬ 
lattice and for that reason all coincide. The value Pxy(O) 
may be decomposed into its intrinsic or coherent and 
extrinsic or incoherent contribution [29, 39]. The later 
part represents in particular all so-called skew scattering 
contributions. Keeping the spin configuration fixed to 
that for T = 0 K (corresponding to a collinear ferromag¬ 
netic ordering) and accounting for the lattice vibrations 
only the resulting p^y{T) rises with temperature T and 
shows a broad maximum around 40 K. For higher tem¬ 
peratures p'^y{T) decreases with T as was found before 
for example in the case of pure Ni in the regime below the 
Curie temperature. [30] On the other hand, if the lattice 
structure is kept undistorted and spin fluctuations are 
accounted for p^y^^{T) monotonously decreases with T 
until the Curie temperature is reached and van¬ 

ishes as the thermally averaged z-component of the mag¬ 
netization {Mz)t vanishes. If both scattering mechanism 
are accounted for one finds again that these do not act in 
a simple additive way. Nevertheless, from the behavior 
of the individual curves p^y(T) and one may 

expect the behavior of the transverse resistivity p^y(T) 
that account for the lattice vibrations and spin fluctua¬ 
tions simultaneously. When comparing these results with 
the corresponding experimental data, a reasonable agree¬ 


ment is found. In particular concerning the presence of a 
maximum for p^y{T) with temperature. The individual 
curves p^y{T) and p^y^^{T) clearly demonstrate the im¬ 
portance of extrinsic contributions to p^y{T) supporting 
the conclusions of Checkelsky et al. [15] concerning their 
role for Fe 2 ,TaS 2 mentioned in the introduction. 


IV. SUMMARY 

Using the KKR-CPA band structure method for dis¬ 
ordered systems the electronic structure as well as mag¬ 
netic and transport properties of ferromagnetic ordered 
and disordered 2H-Fea;TaS2 have been investigated. By 
means of the fully relativistic of the KKR-CPA the spin- 
orbit induced orbital magnetic moment and magneto¬ 
crystalline anisotropy energy could be calculated in par¬ 
ticular. The various magnetic properties were found in 
reasonable agreement with available experimental data 
with clear indications for the strong impact of the degree 
of order on the Fe sub-lattice of the system. In addition, 
a prominent role of correlation effects was found that 
were ascribed to relatively narrow with of the Fe-related 
bands. The temperature-dependence of the longitudinal 
resistivity as well as transverse anomalous Hall resistiv¬ 
ity was studied using the alloy analogy model within the 
framework of the Kubo-Stfeda formalism. The results 
obtained for disordered 2 H-Feo, 28 TaS 2 were found in very 
satisfying agreement with experimental data. Additional 
calculations accounting for thermally induced lattice vi¬ 
brations and spin fluctuations individually clearly showed 
that their contribution to the resistivity is not additive 
and that the temperature-dependence of the longitudinal 
resistivity is nearly exclusively determined by the lattice 
vibrations. 
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